Perturbations of the Kerr spacetime in horizon penetrating coordinates 
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, We derive the Teukolsky equation for perturbations of a Kerr spacetime when the spacetime 

metric is written in either ingoing or outgoing Kerr-Schild form. We also write explicit formulae for 
' setting up the initial data for the Teukolsky equation in the time domain in terms of a three metric 

and an extrinsic curvature. The motivation of this work is to have in place a formalism to study 
the evolution in the "close limit" of two recently proposed solutions to the initial value problem 
in general relativity that are based on Kerr-Schild slicings. A perturbative formalism in horizon 
penetrating coordinates is also very desirable in connection with numerical relativity simulations 
' using black hole "excision". 

<N 
> 

CT) , There is considerable current interest in studying the collision of two black holes, since these events could be primary 
sources of gravitational waves for interferometric gravitational wave detectors currently under construction. On the 
theoretical side, it is expected that the problem of colliding two black holes will be tackled by some combination of 
full numerical and semi-analytical methods. The first three dimensional collisions of black holes are starting to be 
numerically simulated, albeit with very limited resolution and grid size. Long time stability of the codes is also an 
issue. It is therefore of interest to have at hand approximate results which in certain regimes could be used to test 
the codes. Among such approximation methods is the "close limit" approximation [jlj in which the spacetime of a 
I | black hole collision is represented as a single distorted black hole. This approximation has been used successfully to 
£~ I ■ test codes for the evolution of black hole space-times that are axisymmetric and tests are under way for "grazing" 
inspiralling collisions Q . Another realm of application of perturbative calculations is to provide "outer boundaries" 
^ i and to extend the reach of Cauchy codes into the radiation zone far away from the black holes, as was demonstrated 
in ||. Finally, perturbative codes can be used after a full non- linear binary black hole code has coalesced the holes 
to continue the evolution in a simple and efficient fashion as was demonstrated in . 

The perturbative approach requires specifying both a background metric and a coordinate system when performing 
calculations. For evolutions in the time domain such as the ones we are considering, one also has to specify an initial 
slice of the spacetime. In all perturbative evolutions performed up to now the background spacetime has either been 
the Schwarzschild solution in ordinary coordinates or the Kerr spacetime in Boyer-Lindquist coordinates. These 
backgrounds are adequate for instance, for the study of the evolution of "close limits" of the Bowen-York and 
"puncture" |?J families of initial data, which reduce in the "close limit" to those background space-times. 

There have been two recent proposals for alternative families of initial data that have some appealing features 
|^],[)| . Both these proposals are based on the use of the Kerr-Schild form of the Schwarzschild (or Kerr) solutions to 
represent each of the black holes in the collision. Some of these solutions do not have an obvious close limit in which 
they yield a single black hole, although the close limit can be arranged with a simple modification of the original 
proposal. In the case in which the close limit exists, the initial data appear as a perturbation of the Schwarzschild or 
Kerr spacetime, but in Kerr-Schild coordinates. If one wishes to evolve perturbatively the spacetime, this requires 
having the perturbative formalism set up on a background spacetime in Kerr-Schild coordinates. To our knowledge, 
this has never been done in the past. 

The use of Kerr-Schild coordinates appears quite desirable in the context of numerical evolutions of black holes. 
The coordinates penetrate the horizon of the holes without steep gradients in the metric components. This makes 
them amenable to the numerical technique of singularity excision, which is sometimes viewed as the key to long term 
binary black hole evolutions such as the ones needed for gravitational wave data analysis purposes [Tc| ] . 
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Since the Kerr-Schild coordinates are horizon-penetrating, developing a perturbative formalism in these coordinates 
could, in principle, allow the study of perturbations arbitrarily close to the horizon and even inside the horizon. The 
traditional perturbative formalisms, based on Schwarzschild and Boyer-Lindquist (in the Kerr case) coordinates 
cannot achieve this. Having a perturbative formalism that works close to the horizon is desirable in the context of 
the recently introduced "isolated horizons" . One could exhibit perturbatively the validity of several new results 
that are emerging in such context. 

In this paper we will develop perturbative equations in Kerr-Schild coordinates, taking advantage of the fact that 
the Teukolsky formalism is coordinate invariant. We will end by constructing a perturbative equation that is well 
behaved inside, across and outside the horizon. 

The organization of this paper is as follows. In section II we review the Teukolsky formalism (including the setup of 
initial data for Cauchy evolution). In section III we derive the Teukolsky equation in horizon penetrating Eddington- 
Finkelstein coordinates and display a numerical evolution of it. In an appendix we discuss the equation in outgoing 
Eddington-Finkclstein coordinates . 



II. THE TEUKOLSKY EQUATION 

Fortunately, the perturbative formalism due to Teukolsky is amenable to a reasonably straightforward (but 
computationally non-trivial) change of background. The Teukolsky formalism is based on the observation that the 
Einstein equations written in the Newman-Penrose formalism naturally decouple in such a way that one obtains an 
equation for the perturbative portion of the Weyl spinor. In the notation and conventions of Teukolsky (which in 
turn follows those of Newman and Penrose), the resulting equation (in vacuum) is, 

[(A + 3 7 - 7* + 4/i + (i*) (D + Ae-p)- (6* - t* + (3* + 3a + 4tt) (5 - t + 4/3) - 3tp 2 ] ^4 = 0, (1) 

where the quantities in brackets are computed using the background geometry and tp4 is a first order quantity in 
perturbation theory. A similar equation can be derived for the ipo component of the Weyl tensor. We will not 
concentrate on this equation, however, since it has not proven as useful for evolutions in the time domain fl3f| . 
If we now particularize to a background spacetime given by the Kerr metric in Boyer-Lindquist coordinates, 

ds 2 = (1 - 2Mr/E) dt 2 + (AMar sin 2 0/E) dtdip - f dr 2 - £d0 2 - sin 2 9 (r 2 + a 2 + 2Ma 2 r sin 2 0/E) dip 2 , (2) 

where £ = r 2 + a 2 cos 2 and A = r 2 — 2Mr + a 2 (and should not be confused with the Newman-Penrose quantity 
A = n^dfj) and considers the Kinnersley tetrad, 

P= [(r 2 + a 2 )/A,l,0,a/A], (3) 
n" = [r 2 + a 2 ,-A,0,a]/(2E), (4) 
= [iasin0,O,l,i/sin0]/(v / 2(r + iacos0)), (5) 

for which the background Newman-Penrose quantities are, 



p = — l/(r — ia cos 0), 


(6) 


(3 = -p* cot0/(2%/2), 


(7) 


7r = iap 2 sin 0/v / 2, 


(8) 


t = —iapp* sin0/v2, 


(9) 


fi = /V A/2, 


(10) 


j = H + pp*(r-M)/2, 


(11) 


a = n — (3* , 


(12) 


i> 2 = Mp 3 . 


(13) 



and the differential operators D — l^d^, A = fi/'c^, 5 = m^d^, then the resulting equation is the well known Teukolsky 
equation, 
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(14) 



where tfj — (r — ia cos 6) 4, ifj 4. As was discussed in references [jHj, one can write the initial data for the above equation 
in terms of the perturbative three metric and extrinsic curvature. The formulae for the Teukolsky function and its 
time derivative are, 



ip4 = - [Rijki + 2K i[k Kq j ] (1) n i m?n k m l + 8 



(i) 



-4 [Rji - K jp Kf + KKji) {1) n [0 m j] n [0 m l] , 
d t ip4 = ^(0)^0 (^4) - n l m J n k m l [d Rijkl]^ 



(15) 
(16) 



(i) 



+8n [0 m : > ] n k m l 

-An [0 m j] n [0 m l] d Q Rji - 2K p {l d K j)p - 2N (0) K JP K P K^ 
+K jl d K + KdoK jt ] (1) 

+2{ip i (l i A - mi ~5)N l (°) + Mm* - ffH&iN* (0) }, 



where 



ipa = - [Rijkl + 2K i[k K l]:j ] (1) Vn 3 m k n l + 4 



K 



(i; 



(17) 



2 [% - iCfc,.Kf + {l l °n j] m°n l - l [0 n j] n°m l ), 



N(p) = (ffkerr) 1S t ne zer °th order lapse, n l in these equations should be taken to be related to that of the original 
tetrad as n° = N^n®^, n % = n l olig + N l ^n . Latin indices run from 1 to 3, and the brackets are computed to only 
first order (zeroth order excluded). The derivatives involved in the above expressions can be computed in terms of 
the initial data on the Cauchy hypersurface as, 



doK = N (0) K pq K™ ~V 2 N (oh 
d a R = 2K™d K pq + 4N (()) K pq KPK s « - 2Kd K, 

doRtjki = -4iV(o) j/\ ,7 />'/ , - Kj[ k Rqi - -R ( A , ; //; , - Kj^.gi]i) j 

+2gi[ k d Ri]j - 2g j[k d Q R l]i - g l [ k gi ]j d () R + 2K i[k d K l]} - 2K j[k d Q K l]l 



(18) 
(19) 

(20) 



and, 



d K i:j = N, 



(o) 



Rij + KKij - 2K ip KP j - N^ViVjN {0) 



(i) 



(21) 



Remarkably, the above formulae are coordinate independent! Therefore the only adjustment needed to specify 
initial data for the evolution equations we will derive in the next two sections is to insert the appropriate background 
quantities in the above formulae. 



III. THE TEUKOLSKY EQUATION IN KERR-SCHILD COORDINATES: 
INGOING EDDINGTON FINKELSTEIN COORDINATES 

The initial data proposed by [|||] is constructed using ingoing Eddington-Finkelstein (IEF) coordinates (strictly 
speaking, since the initial data might include net angular momentum, one is really talking about the generalization 
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of IEF coordinates to the rotating case, commonly referred to as Kerr coordinates). This is in part due to the fact 
that these families of initial data are currently being evolved using a numerical code where the black holes are treated 
using the "excision" technique. This technique requires coordinates that penetrate the horizon, such as the IEF ones. 
The IEF coordinates (V, r, 0, <p) for the Kerr metric are defined through a redefinition of the time coordinate of the 
Boyer-Lindquist coordinates as, 

V = t + r* (22) 
£ = ip + J |-dr (23) 

where r* is the natural generalization to the Kerr case of the usual Schwarzschild "tortoise" coordinate, and is defined 
by, 

-dr. (24) 



2Mr- 

The codes currently being used to evolve the initial data in IEF are written in terms of a coordinate t — V ~ r. We 
will therefore derive the Teukolsky equation in the (t, r, 0, <p) coordinates. The Kerr metric in these coordinates, 

ds 2 = (1 - 2Mr/£) dP - (1 + 2Mr/£) dr 2 - £d0 2 - sin 2 9 (r 2 + a 2 + 2Ma 2 rsin 2 0/£) dip 2 

- (4Mr/£) dtdr + (4Mra sin 2 0/£) dtd(p + 2a sin 2 0(1 + 2Mr/E) drdtp, (25) 

In addition to changing coordinates, it is immediate to see that one needs also to change tetrads. The usual Kinnersley 
tetrad is singular at the horizon, and therefore leads to a Teukolsky equation that is singular. However, it is easy to 
fix this problem by just rescaling Z M by a factor of A and dividing n M by A. This does not change the orthogonality 
properties of the tetrad, but makes it well defined Q Therefore we re-derive the Teukolsky equation using as new 
tetrad vectors, 

^ = [A + AMr, A, 0,2a] (26) 
^ = 1^,-^0,0] (27) 

This redefinition of the tetrad vectors changes the values of the Newman-Penrose scalars. The new values are, 

e = r — M (28) 
1 r + ia cos 9 

p = — (r + ia cos 9)— (30) 

2-1 

with a, (3, 7r, t, ip2 remain unchanged. The resulting Teukolsky equation reads, 

1 A + 4Alr <9 2 ?/>4 1A<9 2 V>4 ^Mr d 2 ^ 4 a d 2 ijj 4 1 1 d 2 ip A 1 1 d 2 ijj A 

2 £ W ~ 2 £ dr 2 ~£~ drdt ~ £ drdtp ~ 2 £ <90 2 ~ 2 sin 2 0£ dip 2 
3a 2 (r - M) cos 2 9 + r(7r 2 + 4a 2 - HMr) + 4iacos0A \ dip 4 

dr 



r 2 (2r + 11M) + a 2 (2r - 3M) cos 2 - 2m(r 2 + 7Mr + a 2 cos 2 0) cos0 1 dtp. 



:, 4 



£(r — za cos 0) 2 J 9< 

(r 2 + 8a 2 - 9a 2 cos 2 0) cos + 2iar(4 - 5 cos 2 0) 1 di\) A 



2£(r- iacos0) 2 sin0 J 90 

4ar(sin 2 - cos 2 0) - 2i(r 2 + 2a 2 - 3a 2 cos 2 0) cos 1 dip 4 



£(r - iacos0) 2 sin 2 J 9^ 

f 24Mr sin 2 - 19r 2 + (21r 2 - 9a 2 + 7a 2 cos 2 0) cos 2 - 2ia[(7r - 6M) cos 2 - (5r - 6M)] cos 



iacos0) 2 sin 2 I 



^4 = 0. (31) 



1 We wish to thank Steve Fairhurst and Badri Krishnan for suggesting this rescaling. 



This equation can be used to study a variety of different things, including boundary conditions at the excision 
region (inside the horizon), perturbations inside and near the horizon, and most significantly, we could use this 
implementation to compare and perhaps even continue evolutions from full numerical codes that use Kerr-Schild 
coordinates. 

The Penetrating Teukolsky Code (PTC) evolves the following equation, where ip = (r—ia cos(6*)) 4 -04 is the Teukolsky 
function: 



sin 9 d9 V"" " d9 J sin 2 6 d(p 2 ' dtdr drdfi 

4icot6>\ chp 
sin 9 J dip 

- (4r + Ma cos 9 + 6M ) + 2(3 cot 2 9 - esc 2 9)ip = 0. 
at 

To implement this equation numerically, we break the above equation down into a 2+1 dimensional one, using a 
decomposition of the Teukolsky function into angular modes, ip = ^me lra ^ The Teukolsky equation for each m 
mode now looks like the following: 

( S + 2 Mr)^%- A^^-(2aim + 6r-6M)^ (33) 
at 2 or z or 

sm 9 o9 \ 89 J dtdr dt 

+ (4 cot 2 6» - 2 + to 2 esc 2 9 - 4m cot 6» csc $ )?/>,„ = 0. 

We use Lax-Wendroff technique to numerically implement these simplified set of equations exactly as done in []l3f . A 
full discussion of the applications of this code will be presented elsewhere. Here we just highlight in the figures how 
the code indeed evolves perturbations inside and outside the horizon as expected. 



IV. CONCLUSIONS 



We have set up the black hole perturbation framework in Kerr-Schild type coordinates. This framework will be 
useful for comparisons with fully nonlinear numerical codes currently being implemented that run naturally in Kerr- 
Schild coordinates. We have also discussed the numerical implementation of the perturbative evolution equation 
and how to set up initial data in terms of the initial value data that will be available for black hole collisions. 
Implementation of this framework for numerical computations is essentially complete. 

This formalism allows us to study in a natural way perturbations close to the horizon and may also be of interest 
to study and test in a concrete fashion several attractive properties of "isolated horizons" |ll]] . This formalism allows 
us to make several predictions about quantities defined with notions intrinsic to the black hole (like a concept of local 
mass and angular momentum), and their evolution. Several attractive formulae, for instance relating the ADM mass 
to the "local horizon mass" and the "radiation content" can be worked out. Having a perturbative formalism that 
operates correctly near and on the horizon will allow us to test the validity of these formulae. This issue is currently 
under study. 
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Appendix:The Teukolsky equation in Outgoing Eddington Finkelstein coordinates: 



The outgoing Eddington-Finkelstein coordinate form of the Kerr metric might be useful to evolve in a perturbative 
fashion the outgoing wave zone exterior, say, to a tube in which a Cauchy evolution code is used, as in the spirit of 
||. The outgoing Eddington-Finkelstein coordinates for the Kerr metric are obtained by introducing a time variable, 

U = t-r* (34) 
<p = <p - J ^-dr (35) 

where r* is defined as in (p4[). 

As in the previous section, one presumably wishes to write a Cauchy evolution code for the resulting perturbative 
equation. It is therefore appropriate to introduce a time coordinate, 

t = U + r, (36) 

and consider the Teukolsky equation in (t, r, 9, tp) coordinates. The Kerr metric in these coordinates is the same as 
( p5| ) except an opposite sign in the drdt term. 

The Newman-Penrose scalars do not change. The differential operators do. The Kinnersley tetrad in the new 
coordinates is, 

J" = [1,1,0,0] (37) 

r A / 4Mr\ A a, 
n " = ^{ 1 + —)>-^°>^ 

m" = [w sin 0,0, 1, — ^ }/(V2(r + ia cos 9)), (39) 
suit' 

The resulting Teukolsky equation in these coordinates is given by, 

, d 2 ib , d 2 ib d 2 ib 1 d 2 ib d 2 ib d 2 ib 

(E + 2Mr) -i - A-| - -| - -t-^-^I + 4Mr-^ + 2a- 



8i 2 Or 2 dO 2 sin 2 9 dtp 2 dtdr drdtp 

n dib , , dib dib cos 6 dib , \ , ^ 

-2 (2r + M + 2ia cos 9) -L + 2 (r - M) -- - cot 0-±- + 4d — 5--^ + (2 + 4 cot 2 9)ip = 0, (40) 
dt dr d9 sin 2 6 dp v ' Y ' v ' 

which can be rewritten in a form more reminiscent of the ordinary Teukolsky equation as, 

1 d ( . dib\ 1 d 2 ij d 2 t/j d 2 i/j 

smff— =— — — — 2a—— 4Mr— 



sin 9 80 \ 89 J sin 2 9 dip 2 dpdr dtdr 

, dib cos 9 dib , , , 

-2 (2r + M + 2ia cos 9) + Ai— -^-ttt +2(1 + 2 cot 2 6)ib = 0, 



where again ip = (r — ia cos 9) i ip4 : , and we see that the main differences with the equation we derived in the previous 
section are: a) the sign of the terms with the mixed r, <p and r, t derivatives, b) the first order derivative in t and tp 
terms. 

As in the previous cases, the generic initial data formulae we presented are still valid in this case. 
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FIG. 1. Numerical evolution of the horizon penetrating Teukolsky equation. The horizon is at r = 200. We give as initial 
data a Gaussian pulse outside the horizon. The background black hole has a mass of 100 units, and zero spin. Shown above 
are snapshots of the evolution of the m = mode at t = 30, 40, and 50 units. The pulse splits into two pulses, one infalling 
towards the singularity at r = 0, moving smoothly through the horizon. The other pulse moves rightwards and eventually 
escapes to scri. Notice that the speed of light does not appear as a t = r motion in these coordinates! 
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0.8 Initial Data 
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